New Physics in B —B mixing in the light of recent LHCb data 
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We perform model-independent statistical analyses of three scenarios accommodating New Physics 
(NP) in AF — 2 flavour-changing neutral current amplitudes. In a scenario in which NP in Bd—Bd 
mixing and B 3 — B s mixing is uncorrelated, we find the parameter point representing the Standard- 
Model disfavoured by 2.4 standard deviations. However, recent LHCb data on B s neutral-meson 
mixing forbid a good accommodation of the D0 data on the semileptonic CP asymmetry Asl- We 
introduce a fourth scenario with NP in both M^ and T^ , which can accommodate all data. We 
discuss the viability of this possibility and emphasise the importance of separate measurements of 
the CP asymmetries in semileptonic Bd and B 3 decays. All results have been obtained with the 
CKMfitter analysis package, featuring the frequentist statistical approach and using Rfit to handle 
theoretical uncertainties. 
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Flavour physics looks back to a quarter-century of pre- 
cision studies at the B-factories with a parallel theoretical 
effort addressing the Standard Model (SM) predictions 
for the measured quantities jl|. With the parameters 
of the Cabibbo-Kobayashi-Maskawa (CKM) matrix [2| 
overconstrained by many measurements one can predict 
yet unmeasured quantities [3(. Still, the global fit to the 
CKM unitarity triangle reveals some discrepancies with 
the SM, driven by a conflict between B(B — > tv) and 
sin(2/3) measured from Bd — > J/^K 4, 5|. Furthermore, 
in May 2010 the D0 experiment reported a deviation of 
the semileptonic CP asymmetry (dimuon asymmetry) in 
Bd.s decays from its SM prediction [6|,LZ| by 3.2 a |8[. In 
June 2011 this discrepancy has increased to 3.9a [13j. 
In summer 2010 the data could be interpreted in well- 
motivated scenarios with New Physics (NP) in B — B 
mixing amplitudes [4J. In this letter we present novel 
analyses which include the new data of 2011, in particu- 
lar from the LHCb experiment. 

Bq—B q (q = d,s) oscillations involve the off-diagonal 
elements M q 2 and T\ 2 of the 2x2 mass and decay matri- 
ces, respectively. One can fix the three physical quanti- 
ties |M* 2 |, |r* 2 | and 4> q = arg(— Mf 2 /rf 2 ) from the mass 

among the eigenstates, their 



width difference Ar,j ~ 2 |r' 2 | cos<j>q and the semilep- 
tonic CP asymmetry 
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Mf 2 is especially sensitive to NP. Therefore the two com- 
plex parameters A s and A^, defined as 



Mf 2 = M^' q -A 



11 ■■ 



A g ^|A 9 | 



q = d, s, (2) 



difference AA/„ 



2\M? 2 \ 



can differ substantially from the SM value A s = A^ = 
1. Importantly, the NP phases <pf s do not only affect 

a S L, but also shift the CP phases extracted from the 
mixing- induced CP asymmetries in Bd — > J/^K and 
B s -)• J/*0 to 2/3 + <p£ and 2/3 s - (f>f, respectively. In 
summer 2010 the CDF and D0 analyses of B s -*■ J/^cf) 
pointed towards a large negative value of 4>f, while si- 
multaneously being consistent with the SM due to large 
errors. With a large cj>f < we could accommodate DCs 
large negative value for the semileptonic CP asymmetry 
reading Asl = 0.6ag L + 0.4ag L in terms of the individual 
semileptonic CP asymmetries in the Bd and B s systems. 
Moreover, the discrepancy between B(B — > tv) and the 
mixing- induced CP asymmetry in Bd —* J/^>K can be 



removed with (f>^ < 0. The allowed range for <ftf implies 
a contribution to Asl with the right (i.e. negative) sign. 
In our 2010 analysis in Rcf. |4( we have determined the 
preferred ranges for A s and A^ in a simultaneous fit to 
the CKM parameters in three generic scenarios in which 
NP is confined to AF = 2 flavour-changing neutral cur- 
rents. In our Scenario I we have treated A s , Ad (and 
three more parameters related to K — K mixing) inde- 
pendently, corresponding to NP with arbitrary flavour 
structure. Scenario II implements minimal-flavour viola- 
tion (MFV) with small bottom Yukawa coupling entailing 
real A s = A^. Scenario III covers MFV models in which 
A s = Arf is allowed to be complex. In Ref. |4[ we have 
found an excellent fit in Sc. I (and a good fit in Sc. 3) 
with all discrepancies relieved through A d _ s ^ 1, while 
the fit has returned K — K mixing essentially SM-like. 

The recent LHCb measurement of the CP phase $f ^ 
from A™p(B s — > J/^(f>) does not permit large deviations 
of <f>g from zero anymore. This trend was also confirmed 
by the latest CDF results [10|. The current situation 
with the phase 2(j)f^ = — 2f3 s + 4>f and Asl is as follows 
(at 68% CL): 



2<$* = (-32±2i)° D0 

-60° < 2<j>^ < -2.3° CDF [lC 

2<pf* = (-0.1 ± 5.8 ± 1.5)° LHCb J/i^cj) Q 

2(/)f fo = (-25.2 ± 25.2 ± 1.2)° LHCb J/tpf Q 

^sl = (-7.87 ±1.72 ±0.93)- 10" 3 D0 (3) 

Here 2/3 s = 2 sxg(-V ts V?J (V cs V* b )) ~ 2.2° Q. 

From this discussion, there is a conflict between LHCb 
data on B s — > J/ip<j> and the D0 measurement of A$l 
which we cannot fully resolve in our Scenarios I, II and 
III. We therefore discuss a fourth scenario which also per- 
mits NP in the decay matrices rf 2 or Tf 2 . 



RESULTS FOR SCENARIOS I, II AND III 

In Tab. [J we summarise the changes in the inputs 
compared to Tabs. 1-7 of Rcf. 0. Following Ref. 
we have included K^, Ki2, i^n (and the related r de- 
cays) for |V^d| and \V US \- Concerning the measurements 
of (<fi s ,r s ) from B s — > J/ip(j), we have combined the 
CDF and LHCb results by taking the product of their 
2D profile-likelihoods [lfj, LLl| ■ Unfortunately, we could 
not obtain the corresponding likelihood from D0. The 
impact of this omission is mild due to the smaller uncer- 
tainties of the CDF and LHCb results. We have nei- 
ther used the LHCb result on B s — > J/tpfo as only 
4> s (not the 2D likelihood) was provided in Ref. [12J |. 
But we have included the flavour-specific B s lifetime 
t§ s [15j providing an independent constraint on Ar s . 
We analyse the D0 measurement of Agh with the pro- 
duction fractions at 1.8-2 TeV according to Ref. [15| : 



Observable 


Value and uncertainties 


Ref. 


B{K -*■ ev e ) 


(1.584 ±0.020) x 10" 5 


[16] 


B{K -> j«i/ M ) 


0.6347 ± 0.0018 


[12 


B(t -> Kv T ) 


0.00696 ± 0.00023 


[12] 


B(K^ii^)/B(K^n^) 


1.3344 ± 0.0041 


[12 


B(t -»■ Kp t )/B(t ->■ -KV T ) 


(6.53 ±0.11) -10 -2 


[18] 


a 


88.7±t|° 


[14] 


7 


(66 ±12)° 


[14] 


Am d 


0.507 ± 0.004PS" 1 


[16] 


Arris 


17.731 ± 0.045PS" 1 


[22,23] 


A SL 


(-74 ±19) x 10~ 4 


[13] 


0** vs. Ar s 


see text 


[12, n] 




Theoretical Parameter 


Value and uncertainties 


Ref. 






fB. 


229 ± 2 ± 6 MeV 


[14] 






Ib s I!b a 


1.218 ±0.008 ±0.033 


[14] 






Bb s 


1.291 ±0.025 ±0.035 


[14] 






BbJBb, 


1.024 ±0.013 ±0.015 


[14] 






Bk 


(0.733 ± 0.003 ± 0.036) 


[14] 






Ik 


156.3 ± 0.3 ± 1.9MeV 


[14] 






Ik/U 


1.1985 ±0.0013 ±0.0095 


[14] 






a 3 {M z ) 


0.1184±0±0.0007 


[16] 





TABLE I. Experimental and theoretical inputs inputs added 
or modified compared to Ref. [J] and used in our fits. 



f s = 0.111 ±0.014 and f d = 0.339 ±0.031, corresponding 
to A SL = (0.532 ± 0.039)a^ L ± (0.468 ± 0.039)a| L . 

We summarise our results in Tabs. Ql] an d Mil and in 
Fig. Q] (Sc. I) as well as Fig. [2] (Sc. III). Even in Sc. I 
our fit to the data is significantly worse than in 2010 [4J: 
While (j)^ < alleviates the discrepancy of A$l with the 
SM, the LHCb result on (pf^ prevents larger contribu- 
tions from the B s system to Ag^. In Sc. I, we find pull 
values for A$l and <fif — 2(3 s of 3.0 a and 2.7 a respectively 
(compared to 1.2 a and 0.5 a in Ref. [4|). We do not quote 
pull values for Am^s in Sc. I, as these observables are not 
constrained once their experimental measurement is re- 
moved. In contrast to earlier analyses, only one solution 
for A s survives thanks to the recent LHCb determina- 
tion of Ar s > [Hi entailing Re A s > 0. Tab. [TV] lists 
the p-values for various SM hypotheses within our NP 
Scenarios (more information can be found in Rcf. J14|L 



NEW PHYSICS IN Tf 2 OR T? 2 

Several authors have discussed the possibility of a siz- 
able new CP-violating contribution to rf 2 to explain the 
D0 measurement of A$l PUl by postulating new B s de- 
cay channels with large branching fraction. In such mod- 
els also the width difference Ar s typically deviates from 
the SM prediction in Ref. @,0, Ef. Tf 2 is dominated by 



Quantity 


1(7 


3(7 


Re(A d ) 


fl qoo+0.143 
U - OZO -0.095 


0.82l»j 4 


Im(A d ) 


-o.i99ig:8S 


-0.2018:8 


|Ad| 


0.86±g;^ 


o.86l8:l 


<t>* [deg] 


-13.41?:? 


-13.41ft 1 


Re(A s ) 


0.965i°; 3 7 ! 


o.97l8: 3 ° 


Im(A s ) 


-n nn +0 10 

u - uu -0.10 


u - uu -0.32 


|A S | 


fl 077+0121 
U ' J " ' -0.090 


0.9818:? 9 


<£ [deg] 




-oiil: 


cf>f + 2/3 [deg] (!) 


17 +12. 
i '-13. 


iTiS: 


4>f - 2ft [deg] (!) 


-56.8l* 9 


c 7 +66. 


A SL [10- 4 ](!) 


-i5.6i|:2 


-16+ 19 

lu -12 


A S L [lO" 4 ] 


-i7.rl 3 3 : 9 8 


-1812 


a§ L -a| L [lO" 4 ] 


33.61™ 


34±§| 


«sl [lO" 4 ] (!) 


-33.211:? 


qq+25 

°°-13 


«sl [lO" 4 ] (!) 


o.4±S:2 


oi 2 2 ; 


Ar^ps- 1 ] 


0.004801°:°°°™ 


o.oo48i8:88?? 


Ar.fps- 1 ] (!) 


A lcc+0.020 

u.iao_ 079 


n lcc+0.036 
u - iJO -0 .098 


AI^ps- 1 ] 


fl in4+ - 017 
u.iu^_ 016 


A ln4 +0.052 

u.iu^_ 041 


B->ru [lO" 4 ] (!) 


1 Q 41 +0.064 
i.oy:i_ 232 


1 Qzl+0.20 


B -^tv [10^ 4 ] 


1 qc/l+0.063 
1-OJ1-0.095 


i qc+0.19 



TABLE II. CL intervals for the results of the fits in Scenario I. 
The notation (!) means that the fit output represents the indi- 
rect constraint with the corresponding direct input removed. 



Quantity 



Deviation wrt 
SM Sc. I Sc. II Sc. Ill 



tf + 20 

t>t - 2ft 



2.7 a 2.1(7 2.7 a 1.2 a 
0.3 a 2.7 a 0.3 a 2.4 a 



Am d 
Am 3 

Asl 
asL 

Ar s 



0.0 a 


- 


0.0 a 


- 


1.0 a 


- 


1.0 a 


0.9 a 


0.0 a 


- 


1.0 a 


1.3 a 


3.7 a 


3.0 a 


3.7 a 


3.0 a 


0.9 a 


0.3 a 


0.8 a 


0.4 a 


0.2 a 


0.2 a 


0.2 a 


0.0 a 


0.0 a 


0.4 a 


0.0 a 


1.0 a 



B(B -> tv) 



2.8 a 1.1a 2.8 a 1.7 a 



B(B -> tv), Asl 
0f - 2ft, A SL 
S(B -> rt/), <pf - 



4.3 a 2.8 a 4.2 a 3.4 a 

3.3 a 2.7 a 3.3 a 3.2 a 

2ft, A SL 4.0 a 2.4 a 3.9 a 3.2 a 



TABLE III. Pull values for selected parameters and observ- 
ables in SM and Scenarios I, II, III, in terms of the number of 
equivalent standard deviations between the direct measure- 
ment and the full indirect fit predictions. 
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FIG. 1. Complex parameters A d (up) and A s (down) in Sce- 
nario I. Here a cxp = a — 4> d /2. The coloured areas represent 
regions with CL < 68.3 % for the individual constraints. The 
red area shows the region with CL < 68.3 % for the combined 
fit, with the two additional contours delimiting the regions 
with CL < 95.45 % and CL < 99.73 %. The p-value for the 
2D SM hypothesis A d = 1 (A s = 1) is 3.0 a (0.0 a). 



the CKM- favoured tree- level decay b — >• ccs. Any com- 
petitive new decay mode will increase the total B s width, 
which LHCb finds as T s = 0.657 ±0.009 ±0.008 [ll[, im- 
plying r s /r d = 0.998 ± 0.014 ± 0.012 in excellent agree- 
ment with the SM expectation < T s /r d - 1 < 4 • 10~ 4 



211 ] . The new interaction will open new b — >■ s decay 
modes affecting precisely measured inclusive Bd and B + 
quantities |4J. Furthermore, new decays mediated by a 
particle with mass M > Mw will add a term of or- 



< 
E ° 




o 1 

Re A 



FIG. 2. Constraint on the complex parameter A = A d = A s 
from the fit in Scenario III with same conventions as in fig. [l] 
The p- value for the 2D SM hypothesis A = 1 is 2.1 cr. 



Hypothesis 


Sc. I Sc. II Sc. Ill 


ImA d = 
ImA s = 


3.2a 2.6a 
0.0a 


A d = l 

A s = 1 
ImA d = ImA s 


3.0a 0.6a 2.1a 
0.0a 
= 2.8a 


A d = A s = 1 


2.4cr 



TABLE IV. p-values for various Standard Model hypotheses 
in the framework of three NP Scenarios considered. These 
numbers are computed from the x 2 difference with and with- 
out the hypothesis constraint, interpreted with the appropri- 
ate number of degrees of freedom. 



der M^/M A 



^SM.s 



to rf 2 /r°2 1 , while A, 



lw/M- 



normally receives 
In models in- 



a larger contribution of order M, 2 
volving a fermion pair (/, /) in the final state, e.g. those 
with an enhanced B s — > tt decay [19( , one can solve this 
problem through chirality suppression. The extra con- 
tribution to M* 2 is down by another factor of rril/M 2 , 
while that to Tf 2 is affected by the milder factor of 
TO f/ TO b- Quantities like r^ s will not be chirality sup- 
pressed. Therefore it seems not possible to add large NP 
effects to Ff 2 . 

Phcnomcnologically it is thus much easier to postulate 
NP in Tf 2 rather than Tf 2 , because Tf 2 is constituted by 
Cabibbo-suppressed decay modes like b — > tied. Also here 
chirality suppression is welcome to avoid problems with 
Mf 2 , but inclusive decay observables like the semilcp- 
tonic branching fraction or the unmeasured Ar^ pose no 
danger. Clearly, testing this hypothesis calls for a bet- 
ter measurement of ag L . We have studied a Scenario IV 
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value 

1.0 



0.9 
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lm(S d ) 



FIG. 3. Constraints on Im5 d ,Im(5 s in Scenario IV. The ID 



68%CL intervals are lm 3d 



0-92ii;69, lraS s 



1.2 



+ 1.6 



The 



p- value for the 2D SM hypothesis lm5 d - 0.097, lm5 s = 
-0.0057 is 3.2 a. 



d s 

including the possibility of NP in T 12 . We stress that 
Sc. IV permits NP in the |AF| = 1 transitions contribut- 
ing to rf 2 , but not in other |A_F| = 1 quantities entering 
our fits, such as B(B — > tv). Further no new CP phase 
in b — > tics, which would change <fids> ^ s considered. Such 
a phase might further increase the hadronic uncertainty 
from penguin pollution, which is not an issue in the SM 
at the current levels of experimental precision. 
Handy new parameters are 



r 12 /M£ 



ReiT^yMf^y 



Re S q , Im S q amount to 



.s, 



(4) 



(Ar g /AA/ g )/(Ar^ M /AA^ M ) 



and -a!? L /(ArS M /AM g SM ), respectively. The best fit 
values of the SM predictions are 5^ M = 1 + 0.097 i and 
<5f M = 1 — 0.0057 i. ReSd is experimentally only weakly 
constrained. We illustrate the correlation between Im^ 
and Im^ s in Fig. 03 relegating correlations of Re S s with 
lmS dtS to Ref. |l||). The p- value of the 8D SM hypoth- 
esis A d = A s = 1, S d , s = (5™ is 2.6 a. 

We stress that too large values for \S S — <5f M | are in 
conflict with other observables as explained above. We 
have also studied Scenario IV without NP in the B s sec- 
tor (A s = 1 and 5 S = 6 s ,sm)- It could accommodate the 
main anomalies by improving the fit by 3.3er, but with 
large contributions to T^ 2 : Im^ = 1.60^J'°g. 



CONCLUSIONS 

We have performed new global fits to flavour physics 
data in scenarios with generic NP in the B d — B d and 



B s —B a mixing amplitudes, as denned in Ref. |4|. Our 
results represent the status of the end of the year 2011. 
Unlike in summer 2010 the two complex NP parameters 
Ad and A s (parametrising NP in M 12 ) are not suffi- 
cient to absorb all discrepancies with the SM, namely the 
D0 measurement of Asl and the inconsistency between 
B(B -► tv) and A™£{B d ->■ J/^K). Still in Scenario 
I, which fits Ad and A s independently, we find the SM 
point Ad = A s = 1 disfavoured by 2.4 er; this value was 
3.6cr in our 2010 analysis |4j We notice that data still 
allow sizeable NP contributions in both Bd and B s sec- 
tors up to 30-40% at the 3er level. The preference of Sc. I 
over the SM mainly stems from the fact that B(B — >• tv) 
favours 0^ < which alleviates the problem with Asl- 

In order to fully reconcile Asl with cpf^ we have ex- 
tended our study to a Scenario IV, which permits NP 
in both M 12 and T 12 . While this scenario can accom- 
modate all data, it is difficult to find realistic models in 
which the preferred NP contributions to rf 2 (composed 
of Cabibbo-favoured tree-level decays) comply with other 
measurements. There are fewer phenomenological con- 
straints on the Cabibbo-suppressed quantity Tf 2 ; a pos- 
sible conflict with Mf 2 can be circumvented with chirality 
suppression. NP in Mf 2 and Tf 2 with the B s system es- 
sentially SM-likc appears thus as an interesting possibil- 
ity, requiring only a mild statistical upward fluctuation in 
the D0 data on Asl ■ Clearly, independent measurements 
of ag L , a| L and/or <Zg L — ag L are necessary to determine 
whether scenarios with NP in Tf 2 and/or Tf 2 are a viable 
explanation of discrepancies in AF = 2 observables with 
respect to the Standard Model. 
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viding us with the 2D profile likelihood functions needed 
for our analyses. A.L. is supported by DFG through 
a Hciscnbcrg fellowship. U.N. acknowledges support by 
BMBF through grant 05H09VKF. 
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